THU-94-12 



comp-gas/9501001 



Algebraic Spatial Correlations and 
Non-Gibbsian Equilibrium States]^ 

M.H. Ernst and H.J. Bussemaker 



Institute for Theoretical Physics 

University of Utrecht 
Princetonplein 5, P.O.Box 80.006 
3508 TA Utrecht, The Netherlands 



Abstract 

Non-Gibbsian stationary states occur in dissipative non-equilibrium 
systems. They are closely connected with the lack of detailed balance 
and the absence of a fluctuation-dissipation theorem. These states 
exhibit spatial correlations that are long ranged under generic con- 
ditions, even in systems with short range interactions, provided the 
system has slow modes and some degree of spatial anisotropy. In this 
paper we present a theory for static pair correlations in lattice gas 
automata violating detailed balance, and we show that the spatially 
uniform non-Gibbsian equilibrium state exhibits long range correla- 
tions, even in the absence of an external driving field. 

Key words: long rang spatial correlations, violation of detailed bal- 
ance, fluctuation-dissipation theorem, lattice gas automata. 
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1 Introduction 



As highlighted in a recent review by Dorfman et al. a major theme in 
non-equihbrium statistical mechanics during the past few decades has been 
the question under what conditions the correlations in fluids, that consist of 
molecules with short-range interactions, can become long-ranged. The ex- 
istence of generic long range spatial correlations in non-equilibrium station- 
ary states of condensed matter (e.g. systems with a constant temperature 
gradient, driven diffusive systems 0) is by now well understood, and is inti- 
mately connected with the existence of long time tails in Green-Kubo type 
time correlation functions. In this paper we discuss the occurrence of such 
correlations in lattice gas automata (LGA's). 

The long range correlations can be described at three different levels of mi- 
croscopic detail: (i) The macroscopic phenomenological equation (diffusion, 
Navier-Stokes) with a Langevin noise term added Here the transport 

coefficients and noise strengths are phenomenological input in the theory. 
A summary of this approach will be given in section 2. (ii) The kinetic 
(Boltzmann) equation to which a Langevin noise term is added Here the 
transport coefficients are determined by the theory, but the noise strength 
is phenomenological input. This intermediate level of description will not be 
pursued here. For applications to LGA's we refer to Ref. 0. (iii) A fully 
microscopic statistical mechanical description. Here both the noise strength 
and the transport coefficient are calculated from the theory. 

To obtain a description at the fully microscopic level, one starts from the 
equations of motion for the dynamic variables or the phase space density, and 
then derives a BBGKY-hierarchy for the distribution functions. After ap- 
propriate approximations one obtains kinetic equations for the single particle 
and pair correlation functions. Long range spatial and temporal correlations 
are created by sequences of ring collisions. 

If a classical fluid is in thermodynamic equilibrium, long range spatial 
correlations only occur under very special conditions, when parameters like 
the temperature are tuned to a critical value. In general however the corre- 
lation length is flnite. The situation is different when external driving flelds 
or constraints imposed by external reservoirs prevent the fluid from achiev- 
ing a thermodynamic equilibrium state satisflng the constraints of detailed 
balance. In this case generic long range spatial correlations exist in the 
non-equilibrium stationary state, for almost all parameter values. 
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The same concepts are also directly relevant for LGA's that violate the 
constraints of (semi)-detailed balance. Recently, a microscopic theory at the 
level of pair correlations (ring kinetic theory) has been applied to LGA's 
to calculate the equal-time position and velocity correlations in spatially 
uniform equilibrium states 0. The primary goal of the present paper is to 
apply the ring kinetic theory to demonstrate the existence of generic long 
range spatial correlations in uniform equilibria of LGA's without detailed 
balance. This will be worked out in section 3. The coarsest description, by 
means of a fluctuationg diffusion equation, is conceptually the simplest, but 
it captures all essential features. It will be discussed in section 2. 

All three levels of description share the same essential ingredients for the 
existence of such algebraic correlations: (i) the existence of local conservation 
laws or slow (diffusive or hydrodynamic) modes; (ii) lack of detailed balance; 
(iii) some degree of anisotropy, either due to the underlying lattice, or due 
to an external driving field or gradient imposed by external reservoirs. 

2 Langevin Equation 

A clear presentation of the basic theory of non-equilibrium Langevin models 
has been given by Grinstein et al. Q . Here we only indicate the essential 
steps. Consider a locally conserved density h{r,t) = {h{r,t)). On average it 
satisfies a diffusion equation with an anisotropic diffusion tensor. To account 
for the fluctuations, a Langevin noise term is added to this equation. It reads 
after Fourier transformation, 

dtUq^t) = {D^ql + Dyql)h{q,t) +fi{q,t), (2.1) 

where the Gaussian random noise is determined by 

{f,{q,m-q,t')) = 2{B,ql + Byql)5{t-t'). (2.2) 

We allow for spatial anisotropy by having two constants and By. The 
diffusion coefficients and noise strengths B^ are phenomenological co- 
efficients; they are not provided by the theory. Without loss of generality 
we use a two-dimensional presentation to denote longitudinal and transverse 
directions with respect to the wave vector q. 
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The equal time density- density correlation function in the stationary 
state, G{r, oo) = G{r), is given by 



g{r) 



dq 



^tqr 



(2.3) 



Here d is the number of dimensions and the stationary susceptibility can be 
calculated from and ( p.2|) as, 



X{q) = hm {\h{q,t)\ ) - ^ 



i— >oo 



(2.4) 



D{q) D^ql + Dyqf 
where = Qa/Wl is a Cartesian component of a unit vector {a = x,y, . . . ,d) 



Possible scenario's: 



(i) If the system satisfies detailed balance, i.e. if the equilibrium distribu- 
tion is given by the Gibbs distribution, the equilibrium value of the 
susceptibility x{q) — Xo as g — is known from thermodynamics, and 
transport coefficients and noise strengths are related by the fluctuation- 
dissipation theorem, B{q) = XoD{q). The correlation function is short- 
ranged, G{r) ~ xo'^(^)- 

(ii) However without the detailed balance constraint, the equilibrium state 
is not a thermodynamic state, and the corresponding susceptibilities 
are unknown. There is no fluctuation-dissipation theorem imposing 
a relationship between B{q) and D{q). The long wavelength limit 
Xo(^) of the susceptibility is in general anisotropic and depends on the 
direction of q. A rescaling of the integration variables in ( p.3|) , namely 
q = k/r, shows that G{r) ~ r~'^ at large distances. As diffusive modes 
are correlated over time intervals t ^ r"^, the spatial correlations ~ r"*^ 
have an intimate connection with the long time tails in the velocity 
and other current-current correlation functions For later reference 
we quote the explicit form of the correlation function ( p. 3D at large 
distances in two dimensions, where the susceptibility is given by (|2.4| ), 
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As we shall see in section 3, this scenario is realized in LGA's where 
the collision rules do not satisfy detailed balance and break the lat- 
tice symmetry by having different transition probabilities in x- and 
^/-direction. 

(iii) Different scenario's are realized when the phenomenological B{q) and/or 



D{q), and therefore also x{q) i^i (P^) , have a weaker anisotropy as 
q ^ 0, that only shows up in terms of relative order O(g^) or 0{q^). 
Then G{r) decays again algebraically on account of (|2.3| ), with a tail 
proportional to 1/ r'^'^'^ or 1/ r'^"'"^ respectively. 

Suppose that one wants to describe a system where the microscopic transition 
probabilities have the symmetry of a d-dimensional lattice, but where there 
is no externally imposed direction in which the lattice symmetry is broken. 
Examples are the spin-spin correlations on hypercubic lattices as discussed 
by Grinstein [Q, or several non-detailed balance LGA's on square, triangular, 
or FCHC lattices (see sections 3 and 4). 

The general form of the noise strength for small q is 

q'^B{q) = qaqpBap + qaq^q^qsE^f^^g + 0{q^), (2.6) 

and a similar expansion holds for the transport coefficient q^D{q). Greek 
indices denote Cartesian components and summation convention is used for 
repeated indices. The coefficients Sq,/?, Dap and -8^/3755 D'^/s-yS respectively 
tensors of rank 2 and rank 4. On all lattices with 4- or 6-fold symmetry, the 
second rank tensors are isotropic, i.e. = Bdap. 

Fourth rank tensors B'^ and D'^ are only isotropic on the two-dimensional 
triangular lattice and the four- dimensional FCHC lattice, and anisotropic 
on all other lattices in two and three dimensions. In the latter cases the 
susceptibility for g — > behaves as xio) = Xo + (f'X2{,q) and G{r) ~ l/r"'+^. 
In the former cases the susceptibility has the form x(q) = Xo + (fx2 + q'^X'i{(i) 
and G{r) ~ l/r"^"*"^. For later reference we calculate the results for a square 
lattice. There X2{q) = Eg + Ea (q'^ + qfj, and G{r) in (|2.3| ) approaches 

G{r) ~ -GEa/nr^ (r ^ 00). (2.7) 

In the microscopic LGA's violating detailed balance, to be discussed in sec- 
tion 3 and 4, the coefficients B{q) and D{q) can be calculated explicitly 
from the ring kinetic theory, and yield the coefficients of the algebraic tails 
in G{r). 
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3 Ring Equation for LGA's 



In this section we present the ring kinetic theory for velocity and spatial 
correlations in LGA's violating detailed balance with respect to the Gibbs 
distribution. Such correlations have been observed in several models by var- 



ious authors R 10, 11 



Starting from a BBGKY hierarchy for the n-particle distribution func- 
tions, we use cluster expansion techniques to derive approximate kinetic equa- 
tions. In zeroth approximation the standard nonlinear Boltzmann equation 
is obtained; the next approximation yields the ring kinetic equation, similar 
to that for hard sphere systems, describing the time evolution of pair cor- 
relations. The derivation of the generalized Boltzmann equation and ring 
equation can be found in Ref . [0] . Here we only sketch the results in so far as 
they are needed for the present analysis. Similar cluster expansion techniques 
have been developed in Ref. [Q. 

The collision step of a LGA is defined by the strictly local transition 
probabilities Ago- from a precollision state s at a single node to a postcoUision 
state a at the same node. It is followed by a propagation step. 

Much is known about LGA's that satisfy the condition of semi-detailed 
balance, J2s^sa = 1, or (stronger) detailed balance, Asa- = A^s, with respect 
to the Gibbs distribution. The equilibrium Gibbs distribution depends only 
on the globally conserved quantities, such as the total number of particles A'", 
the total momentum P or the total energy E. The Gibbs distribution is com- 
pletely factorized: no equal-time correlations between occupation numbers 
exist. 

In LGA's without detailed balance J2s^sa 7^ 1- Starting from a com- 
pletely factorized precollision state, the collisions will then create postcoUi- 
sion correlations between occupation numbers at the same node. The post- 
collision on-node correlation for i ^ j are given by []^, |TT 



(cTiOTj)* = J2 cr^ajAsaF{s) 7^ fifj, (3.1) 

where we have introduced the average occupation number /j = (sj). The 
completely factorized distribution function F{s) is given by 

Fis) = l[ftil-0-"- (3.2) 



6 



Due to subsequent propagation and collision steps the postcoUision on- 
node correlations will be transformed into on- and off-node precollision cor- 
relations, via a process of scattering by other particles in the system (ring 
collisions) . 

If we neglect all correlations between occupation numbers, the time evo- 
lution of the one-particle distribution function, /j(r, t) = (sj(r,t)), is given 
by the nonlinear Boltzmann equation 0, 

Mr + c„ t + 1) = Mr, t) + t)), (3.3) 

where r labels the nodes and Cj the nearest neighbor links. For the detailed 
definition of fij(/) in terms of Asa we refer to Ref. [0]. At a more refined 
level of description we include all pair correlations between fluctuations {6si = 

G,,{r,r',t) = {6s,{r,t)6s^{r',t)). (3.4) 

The coupled time evolution of fi{r,t) and Gij{r,r' ,t) is then given by the 
generalized Boltzmann equation 

Mr + c„ t + 1) = Mr, t) + a(/(r, t)) + ^ ^.^(/(r, t))G'fc,(r, r, t), (3.5) 

k<e 

together with the ring kinetic equation 

Gij{r + d, r' + Cj, t + 1) =6{r, r')Bij{r, t) + ^iiM^kiir, r', t). (3.6) 

k,e 

Here the linear pair collision operator is given by 

^ij,ki = {Sik + ^ik{f{r, t))}{6ji + ^ji{f{r, t))}. (3.7) 

In the equations above we have introduced flik = dQi/dfk and ^li^ke = 
d'^Qi/dfkdfi. The on-node source term Bij{r,t) depends on fj{r,t) and 
Gij{r,r,t), and essentially represents the correlations created by the colli- 
sion step. 

The explicit form of the matrix elements Q are not needed in the present 
analysis, where we focus on the large-r behavior of the solution of (|3.6|) in 
the uniform equilibrium state. The on-node source term contains a dominant 
contribution, referred to as the simple ring approximation, which will be given 
below for the special case of equilibrium. Here we neglect subleading terms. 
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linear in the pair correlation function 0, which correspond to the 'repeated 
ring approximation' in the theory of continuous fluids. 

A spatially homogeneous equilibrium solution obeys the relation fi{r, oo) 
= fi and Gijir, r', oo) = Qij{r — r'). For a given average occupation /j, the 
correlations can be obtained from a linear equation, 

with vq the volume of the unit cell of the direct lattice (fo = 1 for square and 
cubic lattices and Vq = |"\/3 for triangular lattices), and the pair streaming 
operator given by 

Sijiq) = exp[iq ■ (cj - c^)]. (3.9) 
The on-node source term 

Bij = Bij{r, oo) = {(JiCTj)* - {siSj) (3.10) 

equals the change in on-node correlations caused by the collision step in 
equilibrium, where (...)* is defined in (|3.1j ) {Bij = in SDB-models). 

The spatial correlation functions of most interest are those between con- 
served densities, defined as 

Gatir) = Y^a,h,g,j{r) = {ab\g{r)), (3.11) 

where aj and 6, are collisional invariants. In a purely diffusive system Oj = 
bi = 1, in a thermal fluid-type LGA a^, hi = {1, Cai, \c^}- Moreover, we note 
the important conservation law {ab\B) = as implied by ( p.lOj ), stating that 
the source term B is orthogonal to products of collisional invariants. 

The numerical calculation of ^ij (O) for given fi involves matrix inversion 
and integration over the first Brillouin zone. Since fi again depends on Qij{0), 
we use an iterative scheme to find a self-consistent solution {fi, ^^(0)}. Once 
GijiO) is known, all off-node correlations Qij{r) can be calculated from it in 
a straightforward way. For details we refer to Ref. 0. 

The analytical calculation of ( p.8|) can be carried through for large spatial 
distances (r —>■ oo) or small wave numbers {q 0). The essential observation 
is that the pair operator s{q)uj in (|3^ ) has slow (diffusive) modes Xtiv{.Q) 
with eigenvalues behaving like Kp^p{q) ^ 1 — q^D^y{q) + . . . for small g, where 



1 — s[q)uj 



s{q)B 



(3.8) 
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D^y{q) may depend on the direction of q. The slow modes give rise to singular 
denominators in ( p.8|) , and are responsible for the large-r behavior of Qab{f)- 
To analyze this more systematically it is convenient to make a spectral 
decomposition of the pair operator s{q)uj in ( p.8|) , i.e. 



lab{r) ^ VqY^ 



f dg ^ifj,r {ab\x^^u{q)){Xf.uiq)s{q)\B) 



Here x^ll' cind X/iv cire respectively right and left eigenfunctions (product 
modes) of the pair operator s{q)u! with eigenvalue A^^{q). These prod- 
uct modes form a complete bi-orthogonal set of basis functions in the pair 
space, i.e. {XnulXfi'u') = ^ ij.fi' ^uu'- The asterisk in ( |3.12D indicates that the fiu- 
summation is restricted to pairs of slow (hydrodynamic or diffusive) modes. 

To make these results more explicit we specialize to a purely diffusive 
model on a square lattice. Here = 6j = 1 are the only collisional invariants, 
and the only slow pair mode Xoniq) is a product of single particle diffusive 
modes. The conservation laws imply = 0, and the lattice symmetry 

guarantee that {xDD{q)s{q)\B) = q^B{q) where B{q) has the general form 
( [2.61 ). The denominator, 1 — Kooiq) = 2q'^-D(qr), has a similar form for small 
q. Moreover, one can show that {11\xdd) = l + O^q"^). The (non)-isotropy of 
the second and fourth rank tensors D and B^, respectively, depends 
on the microscopic collision rules of the lattice model under consideration. 

In summary, we have established that the spatial correlation functions 
Qab{f) in equilibrium approach for large r the generic form ( |2.3| ) and ( |2.4| ) of 
the phenomenological theory of section 2, and in addition we have provided 
expressions for the tensors B, D and B^, D^, leading to a complete micro- 
scopic theory for long range behavior of the spatial correlation functions of 
LGA's violating the detailed balance constraints. 



4 Interacting Random Walkers 
4.1 Spatial Symmetry 

To illustrate the general results of the previous sections we first consider a 
model of interacting random walkers on a square lattice with Fermi exclusion. 
There is no breaking of the lattice symmetry in the microscopic transition 
probabilities. The model allows for at most one particle per state {r, Cj} 
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{i = 1,2,3,4) where r labels a lattice node and Cj a nearest neighbor link. 
The in-state or pre-transition state of the node is described by the set of 
occupation numbers s(r) = {sj(r);i = 1,2,3,4}, and similarly the out-state 
or post-transition state is described by a{r) = {crj(r);i = 1,2,3,4}. The 
dynamics consists of local transitions, followed by propagation. The local 
transitions conserve the number of particles. 

If there is only one particle at node r, it behaves like a random walker; 
if there are three particles at r, or equivalently one hole, the hole behaves 
like a random walker. If there are no particles or no holes, in- and out-states 
are identical. To describe the two-particle interactions, let {ij} denote the 
two-particle state in which links i and j are occupied. Then the transitions 
are defined as, 

a 

{13}, {24} ^ {12}, {23}, {34}, {41} (4.1) 

where a denotes the transition probability of any state on the left to any on 
the right, and (3 the probability for the reverse transition. We note that the 
states on the right, and those on the left of ( [4.1| ) are not connected by lattice 
symmetries. If we denote the transition probabilities among the four states 
on the right hand side by 7, and those among the ones on the left by 6, then 
we have the normalization conditions 4a + 25 = 1 and 2/5 + 47 = 1. 

In case a ^ (3 the matrix of transition rates, Asa 7^ ^o-s, does not obey 
the (semi)-detailed balance conditions of section 3. The transitions in (|4.1| ) 
create correlations between different velocity links, and the pre-transition 
correlations (sjSj) differ from the post-transition ones {aiaj)* in ( p.lO| ), even 
in the uniform equilibrium state. Consequently, Bij in (|3.10|) is nonvanishing. 



and so is q B{q) in ( |2.6| ). Because of the square symmetry the second rank 
tensor Ba/3 = B6a/3 in (p.6|) is isotropic, but the fourth rank tensor B^j^^g 
is not, which leads immediately to the 1/r^ tail in ( |2.7| ). The coefficient Ea 
is proportional to (a — — /)^ where / = |p is the reduced average 

density. 

In case a = f3 the detailed balance constraints are satisfied, and the tran- 
sitions do not create any correlations in the stationary state, i.e. (aiaj)* = 
{siSj) = fifj, and long range correlations are absent. The details of these 
calculations and the explicit value of the coefficient Ea of the algebraic tail 
will be given elsewhere |jT3[ . 
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4.2 Broken Spatial Symmetry 



The long range correlations will be much stronger if the microscopic transition 
probabilities are anisotropic. It is a driven diffusive system that remains 
spatially uniform in the stationary state. For purpose of illustration we 
consider again a diffusive 4-bit LGA on a square lattice where the local 
transition probabilities in x- and y-directions are different, i.e. 

A,, = -^S{p{s),p{a))exp[bj,{s)j,{cr)]. (4.2) 
Z{s} 

Here jx{s) = J2i CxiSi{r) represent the local particle current in the x-direction. 
The local number of particles p(s) = ^^^^^(r) in the in-state, as well as p{a) 
in the out-state is conserved, and Z{s) is a normalization constant chosen 
such that J2a^scr = 1- For positive 6- values transitions with parallel in- and 
out-currents in the x-direction are favored; those with anti-parallel currents 
are suppressed. There is no bias in the y-direction. If 6 = the LGA 
satisfies detailed balance and represents random walkers interacting through 
the Fermi exclusion rule. All spatial correlations are vanishing. If 6 7^ 
there exists a spatially uniform equilibrium state with long range correlations 
~ The qualitative features of this model are similar to those of the 

diffusive lattice gas, studied in 

It is important to emphasize that the constraints of (semi)-detailed bal- 
ance, i.e. 

J2Fo{s)Asa = Fo{a) (SDB) 

s 

Fo{s)A^ = Fo{cr)A^, {DB) (4.3) 

refer to the factorized Gibbs distribution Fq{s) = C exp[ap{s) + ...], which 
depends on the state variable s only through the conserved densities p{s), 
etc. This implies Fo{s) = FQ{a) and the SDB constraints reduce to the 
familiar form J2s ^su = 1- Similarly the DB constraints in ( [4.3|) reduce to 
A —A 

Therefore the model properties J2s -^sa 7^ 1 or A^^ ^ A^-s do exclude the 
existence of the Gibbsian equilibrium state, but do not necessarily imply the 
existence of position and velocity correlations in equilibrium. This statement 
can be slightly rephrased by saying that the model properties J2s ^sa 7^ 1 or 
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Asa 7^ do not exclude the existence of a completely factorized equilibrium 
distribution. We present an example. 

Consider a variation on the LGA in ( ^.2| ) of a driven diffusive system, 
where 

Asa = -^S{p{s),p{a))exp[j{(T) ■ E]. (4.4) 

Here E is an external bias field and j{s) = J2i CiSi{r) the local (non-conserved) 
particle current. Then clearly, Ylis^sa 7^ 1- Nevertheless, the SDB- or DB- 



condition (4J) does admit a completely factorized (normalized) equilibrium 
solution, 

Fois) = Cexp[«p(s) + jis) ■ E], (4.5) 

which can also be cast in the standard form (|3.2| ) with an equilibrium single 
particle distribution of Fermi-type, 



fi = [l + exp(-a -Ci-E)] \ 



(4.6) 



Establishing the stability of this equilibrium state (|4.5|) - ([4.6|) requires an H 



theorem. In fact, at this meeting H. Chen |Ti| has presented a proof of the 
H-theorem for a lattice gas that satisfies the DB-condition ( [4.3|) with respect 
to a factorized equilibrium distribution ( |4.5| ) with exp[—j{s) ■ E] replaced by 
Ilj{Fj/ Ij), where Fj and Ij are arbitrary positive constants. 

Concluding, we can say that the model satisfies the SDB-constraints ( |4.3| ) 
with respect to the factorized distribution Fq{s). The source term Bij in 
( |3.10D vanishes for this model and consequently all on- and off-node correla- 
tions vanish in the equilibrium state. 



5 Conclusion 

We have studied LGA's with local (on-node) interactions that possess slow 
(diffusive or hydrodynamic) modes and lack (semi)-detailed balance (SDB) 
with respect to the Gibbs distribution. As is well known [Q, ^, for long 
times such systems approach non-thermodynamic spatially uniform states, 
which exhibit long range correlations, if the isotropic spatial symmetry is 
broken by externally applied fields or reservoirs. The surprising result is 
that there exist LGA's that exhibit long-range spatial correlations ~ r"*^"^ 
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in the equilibrium state without any fields or reservoirs that break the spa- 
tial symmetry. These correlations are absent in LGA's satisfying (semi)- 
detailed balance. The property J2s^sa 7^ 1, i-e. violation of the standard 
SDB-constraints, does not necessarily imply the existence of long range cor- 
relation, nor does it exclude the existence of a completely factorized equi- 
librium distribution. This was shown by explicit construction of a diffusive 
LGA model that satisfies the SDB-constraints with respect to a non-Gibbsian 
factorized equilibrium distribution (see ( [4.5|) and ([4.6| )). A similar fiuid-type 



LGA has been presented at this meeting by H. Chen flj]. In continuous fiuids 
long-range correlations and non-detailed balance stationary states can only 
exist in open systems, externally driven or constrained by reservoirs, because 
the thermodynamic (Gibbs) equilibrium state satisfies detailed balance. 
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